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We study the superfluid phase transition in the two-dimensional (2D) excitonic system. Employing 
the extended Falicov-Kimball model (EFKM) and considering the local quantum correlations in the 
system composed of conduction band electrons and valence band holes we demonstrate the existence 
of the excitonic insulator (El) state in the system. We show that at very low temperatures, the 
particle phase stiffness in the pure-2D excitonic system, governed by the non-local cross correlations, 
is responsible for the vortex-antivortex binding phase-field state, known as the Berezinskii-Kosterlitz- 
Thouless (BKT) superfluid state. We demonstrate that the existence of excitonic insulator phase is 
a necessary prerequisite, leading to quasi-long-range order in the 2D excitonic system. 


I. INTRODUCTION 

The interesting quasiparticles, excitons, play the fun¬ 
damental role in the physics related to the recent rev¬ 
elations on the excitonic insulator (El) state, EH1 ex¬ 
citonic Bose-Einstein condensation (BEC) [Ml and the 
excitonic superfluidity, [Ml- In analogy with Cooper 
pair condensate [l 4| , one can naturally expect that the 
electron-hole pairs (excitons), being the bosons with neu¬ 
tral charges, should eventually undergo BEC or the Bose- 
superfluid Berezinskii-Kosterlitz-Thouless (BKT) transi¬ 
tion El EH (in the case of negligible quantum dissipa¬ 
tions) at the very low temperatures. The first one is typ¬ 
ical for the three dimensional (3D) system, where it is 
the dominant phase transition, whereas, the second type 
of transition is typical for two-dimensionalJ_2D) systems, 
where the long-range order is absent [171, ill ■ As in the 
case of the usual Bose gases, there should be a relation 
between these two excitonic phase transitions as it is re¬ 
ported recently EH- It is well known that the BEC in 
the interacting uniform systems occurs only for D > 2. 
However, the absence of BEC does not necessarily imply 
the lack of a superfluid phase transition in D = 2, assum¬ 
ing that the well-defined conditions are satisfied by the 
system EMI. In sharp contrast to the 3D case, inter¬ 
actions in 2D Bose system cannot be treated as a minor 
correction to the BEC picture of the ideal Bose gas and 
they qualitatively change the behavior of the system. In 
fact, the excitonic superfluidity requires the interacting 
Bose gas, El on the contrary, BEC does not. The BEC 
is always superfluid and the existence of the critical Lan¬ 
dau velocity El in a Bose-Einstein condensed gas is well 
known j2Cj. On the other hand, the ideal Bose gas could 
exhibit the phase transition to the BEC state being ab¬ 
solutely not superfluid. 

Despite the continuous attempts, using the ultra-high 
quality materials, to observe superfluidity in bilayer 
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electron-hole systems, such as the double wells in the 
GaAs-AlGaAs heterostructures [2ll - i23j and graphene bi¬ 
layers (in the case of graphene, barriers as thin as 1 nrn), 
the superfluid phase has not been observed yet, except 
of quantum Hall regime in the presence of the exter¬ 
nal magnetic field, where the physics is quite different 
0. Accordingly, it may seem that the electron-hole su¬ 
perfluidity at the vanishing external magnetic field will 
never occur in a solid state system, however, it has been 
shown EH that a double bilayer graphene system, sep¬ 
arated by barrier of thickness 1 nrn, should generate an 
excitonic superfluid at experimentally attainable temper¬ 
atures, and in the case of the absence of the external 
magnetic field. One of the key reasons why the excitonic 
superfluidity is so hard to observe experimentally in 2D 
case is in fact related to the dominant role of quantum 
fluctuations at low dimensions and at low temperatures, 
when the very large zero-point oscillations are present. 
This peculiarity is due to the absence of any real heavy 
particle in the electron-hole (e-h) system. 

In the present paper, we address the role of the parti¬ 
cle phase coupling in the purely 2D in-plane interacting 
excitonic system. We explore the low-temperature quan¬ 
tum collective behavior of the excitons and we extend the 
theoretical works mentioned above, by showing that the 
formation of the excitonic superfluid state is governed by 
the non-local cross correlations between nearest neigh¬ 
bors (n.n.) excitonic pairs in contrast to the formation 
of the El state, where the local on-site correlations are 
important. We derive the BKT transition lines, and we 
discuss the values of the physical parameters entering in 
the system. 


II. THE MODEL 

For the study of the El state and the excitonic super¬ 
fluidity in 2D, we have chosen the two-band extended 
Falicov-Kimball model (EFKM), SI due to its large 
applicability for treatment of the electronic correlations. 
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The Hamiltonian of the EFKM model is given by 


( 3 ) 


u = -t c 22 [ c ( r ) c ( r, )+^- c -] - m '22 — 

<r,r'> r 

-tf 22 [/( r )/( r ') + h -°] + ec 2 e/ + 

<r,r') r 

[ n 2 ( r )“” 2 W] • (!) 

r 

Here /(r) (c(r)) creates an / (c) electron at the lattice 
position r, the summation (r, r') runs over pairs of n.n. 
sites of 2D square lattice, . The density type short hand 
notations are introduced n( r) = n c (r) + n/(r) and h( r) = 
n c (r) — rif( r). Next, t c is the hopping amplitude for c- 
electrons and e c is the corresponding on-site energy level. 
Similarly, tf is the hopping amplitude for /-electrons and 
ef is the on-site energy level for /-orbital. The sign of the 
product t c tf determines the type of semiconductor, for 
t c tf <0 ( t c tf > 0) we have a direct (indirect) band gap 
semiconductor. The on-site (local) Coulomb interaction 
U in the last term of the Hamiltonian in Eq. JT]) plays the 
coupling role between the electrons in the / and c sub¬ 
systems. The chemical potential p is p = p — e, where 
e=(e c + e/)/2. We will use t c = 1 as the unit of energy 
and we fix the band parameter values e c = 0 and e/ = — 1. 
For the /-band hopping amplitude tf we consider the 
values tf = —0.3 and tf = —0.1. Throughout the paper, 
we set k B = 1 and h = 1 and lattice constant a = 1. 


III. THE El STATE DISCUSSION 


The Hamiltonian in Eq. m is containing two separate 
quadratic terms and is suitable for decoupling by func¬ 
tional path integration method [251 . We use imaginary- 
time fermionic path integral techniques, and we intro¬ 
duce the fermionic Grassmann variables /(rr) and c(rr) 
at each site r and for each time r, which varies in the in¬ 
terval 0 < t < /3, where /3 = 1/T with T being the ther¬ 
modynamic temperature. The time-dependent variables 
c(rr) and /(rr) are satisfying the anti-periodic boundary 
conditions x(rr) = —x(r r + /3), where x = f or c. After 
decoupling the last interaction term in the Hamiltonian 
in Eq.(|T| we will have for the grand canonical partition 
function of the system 

^GC = j [DcDc] [VfVf] [VV] [ 2 ^] e -S[c,c,/,/W e ], ( 2 ) 


where the action in the exponential is given by 

V 2 (rr) g 2 (rr) 
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—iV (rr) n (rr) — g (rr) n (rr) 


U 
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The new variables V (rr) and g (rr) in the action are 
the decoupling fields for quadratic terms in the Hamil¬ 
tonian, in Eq. proportional to n 2 (rr) and h 2 (rr) 
respectively. Next, Ss[f,f] and Sg[c,c] are Berry ac¬ 
tions for / and c-electrons and they are defined as follows 
Sb[x,x] = X] r Jcf dra;(rr)i;(rr), where x(r r) = d T x(r r) 
is the time derivative. Next, we will factorize usual elec¬ 
tron operators / and c in terms of new fermionic variables 
/ and c coupled to the unitary charge-carrying U(l) ro¬ 
tor. To this end we write the potential V (rr) as the sum 
of a static and periodic part V (rr) = Vq + V (rr) . Then, 
for the periodic part, we introduce the 17(1) phase field 
kp (rr) via the “Faraday”-type relation j2(| 


V (rr) 


dp (rr) 
dr 


( 4 ) 


For the static part Vo and p(rr)-field, the saddle-point 
evaluations give Vq P ' = iHp- — ip, and g s ' p ' = ^ — 
/ . Here n is the average total particle density 
n = {n c (rr)) + (n/ (rr)) (furthermore, we will fix n = 1, 
corresponding to the case of half-filling 01 ) and n is the 
average of the difference of particle densities n = (h (rr)). 
Then the partition function of the system becomes 


Zgc = j [T>SDc\ [VfVf] [Dip\ e ~s[c,cJJM (5) 


and the total action in Eq.® reduces to 

«S[c, c, /, /, <p] = S eS [ip\ + S B [c, c] + S B [/, /] - 

rP 
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rP 
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Here 


5 effM = 22 


dr 


</ 2 (rr) 2 p . 


U 


- a7 0(rr) “ 


—ic/(rr)n(rr)] 


( 7 ) 


is the phase-only action, which contains fluctuating imag¬ 
inary term ic/(rr)n(rr). The chemical potentials p n and 


Un 


Pn are defined as p n = Hp — p and pa = 2 . 

Next, we perform the local gauge transformation to 
new fermionic variables /(rr) and c(rr). For the elec¬ 
trons of / and c orbitals, the U(l) gauge transformation 
could be written as 


x(tt) 
x(r r) 


= %)■ 


x(r t) 
x(r r) 


( 8 ) 
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where U[}p) is the U(l) transformation matrix U(<p) = 
I ■ cos </?(rr) + ia z • sin tp(rr) with the unit matrix I and 
& z being the Pauli matrix. 



0 1 2 3 4 5 6 7 


U/t 


FIG. 1: The excitonic order parameter A as a function of 
Coulomb interaction parameter U normalized to the hopping 
amplitude t. 


After those transformations we can write the total ac¬ 
tion of the system in the Fourier-space representation 
introducing the vector-space notations and, furthermore, 
we will derive gap equation for the excitonic order pa¬ 
rameter A. The effective phase averaged action of the 
system in the Fourier space takes the following form 


<S e ff 




1 

Jn 


E [ 2 kK),/kK) 

k v n 


X 


x (k, v n ) 


CkOn) 

/k On) 


(9) 


Here Q x (k, u n ) is the inverse of the Green function ma¬ 
trix, given by 


Q *(k ,u n ) 


E&n) -A \ 

-A E£(v n ) J ’ 


( 10 ) 


where single-particle Bogoliubov quasienergies E k (z/„) 
and E^{v n ) are given as E£[v n ) = eg- iv n — tk, E[(u n ) = 
e^ — iUn — tk. Next, fk and fk are band-renormalized hop¬ 
ping amplitudes tk = 2tgs7k and tk = 2tgs7k, where gs 
is the bandwidth renormalization factor 


g B 


; -i[v(rr)-V3(r'r)] 


|r—r' | —d 


(ii) 


and 7k is the 2D lattice dispersion 7 k = cos (k x d x ) + 
cos (kydy), with d a ( a = x,y ), being the components of 
the lattice spacing vector d = r — r' with r and r' nearest 
neighbors site positions. For the simple, square plane 
we have d a = a. The calculation of g b within the self- 
consistent harmonic approximation (SCHA) [23, 28) is 
given in the Section [V] The quasiparticle energies e ? and 
eg are of Hartree-type and they are defined in the theory 
by relation e x = e x — y + Uriy + i (<^(rr)) , where y means 
orbital, opposite to x. 

The El low-temperature phase is characterized by lo¬ 
cal excitonic order parameter A = U (c(rr) f (tt)^ . The 

El state develops from local on-site electron-hole correla¬ 
tions. The expectation value, given in the expression of 
local El order parameter,could be calculated in the frame 
of path integral method [25[ as well as the fermion density 
averages of the respective band levels n x = (i(rr);r(rr)). 
We get a set of coupled SC equations for the El order- 
parameter A, single-particle fermion densities n x and El 
chemical potential y 


[/(*£) + /(£k)]=!, 

k 



/(£+)-/(E k ) 
^k+4A 2 


A = - 


U A 

lv~ 


E 


/(£+)-,/(E k ) 
V^ + 4A2 


( 12 ) 

(13) 

(14) 


Here £k = — tk + £g+Ik —ej is the quasiparticle dispersion 
and the energy parameters E k and E k are defined as 

E k = - f^—tk + eg — ?k + £/ ± sjt ;k + 4A 2 ^j . (15) 

The numerical solution of the system of coupled SC equa¬ 
tions Eas. (fl2l) - (fl4l) is performed for fixed value of the to¬ 
tal particle density n = nj+ng = l and k-summations 
were performed with the (100x100) k-points in the First 
Brillouin Zone (FBZ). The finite-difference approxima¬ 
tion method is used in numerical evaluations, which 
retains the fast convergence of Newton’s method f29t | . 
The accuracy of convergence for numerical solutions is 
achieved with a relative error of order of 10~ 7 . In Fig. [T] 
the numerical results for the local excitonic order param¬ 
eter A for the El state are presented. Different values of 
t are considered. 


IV. U(l) PHASE VARIABLES 

In this Section we integrate out the fermions in the 
partition function in Eq.@ and we obtain the bosonic 
total phase action of the system. We will show how the 
non-local (i.e. between the excitonic pairs on different 
lattice sites) fermionic correlations give the main contri¬ 
bution to the phase stiffness of the ensemble of excitons. 
The partition function in Eq.([2]) could be rewritten as 
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FIG. 2: The parameter J, as a function of the Coulomb in¬ 
teraction parameter U normalized to the hopping amplitude 
t. In the inset, we have presented the shape of the density of 
states p(e) in the case of the 2D square lattice. 


-Z'GC = J 


(16) 


where the effective phase action in the exponential is 


= - In / [: DcDc} 


o-s[c,c,f,fM. After ex- 


VfVf 

panding the logarithm up to second order in the cumu- 
lant series expansion (higher terms are not considered), 
we find for the effective action 

SeffM = 5 0 + (^)s eff [c,c,/,/] - 

( S )s e « [B,cJ,/] “ ^Seff [£,£,/,/] 

(17) 


Here 5 e ff 


c, c, /, / 


is the effective fermionic action 


given in Eq. ©■ We have replaced above (...) g j —»• 
(•••)s f f[ac?/]’ ma ^- e the calculations in a self- 

consistent way. We will examine the four-fermionic terms 
in Eq. urn thereby treating relevant part of non-local 
fermionic correlations. Therefore, the important part of 
the effective phase action is 


SeffM = SoM + 5/M, (18) 


where SoM is the phase-only action after U(l) gauge 
transformation 


5 oM = Y 


r <^ 2 (rr) 2 /i . 


dr 


U 


- iu^ rT) 


(19) 


and [(/j] = — \ («S 2 ) 5 ff ^ e J • After calculating all 

averages in the expression of <Sj[<p] and after not com¬ 
plicated evaluations we rewrite the action 5,/[yj] in the 
form 


Sj W\ 


P 

dr J(rr, r , r / ) x 

<r,r'> 

x cos 2 Yp{yt) — cp(r V)], 



( 20 ) 


where the exciton phase stiffness parameter J is given by 
the relation 


J = 


A 2 tt 


E 

\c' 


V& + 4A 2 


Ai(k, k') tanh 



—A 2 (k, k 7 ) tanh 



( 21 ) 


The parameters Ai(k,k') and A 2 (k, k') in Ea. (l2Tl) are 
defined as Ai(k, k') = (E/ — E /,) 1 ■ (E/ — 1 and 

A 2 (k, k') = ( E ,7 - E^,) 1 • ( E ^ - E£,) '. The form of 
J in Eq. m indicates that the phase stiffness in the sys¬ 
tem of excitonic pairs is characterized by an energy scale 
proportional ~ (A t e th)/(t e + th) for all the values of the 
Coulomb interaction parameter U and it is related to the 
motion of the center of mass of e-li composed quasipar¬ 
ticle, because {t e t h )/(t e +t h ) ft (m e +m h )~\ applying 
that the exchange coupling parameter becomes propor¬ 
tional to the excitonic_BEC-like quasi-condensate critical 
temperature in 2D [30]. In numerical evaluation of the 
stiffness parameter J, given in Eq . (1211) we transform the 
summations over k into energy integrals, by introducing 
density of states p(e) = -k X)k <5 (e — 7 k)- For the tight- 
binding hopping matrix on the square lattice we have 
p(e) =K(1- e 2 /4) / 7 r 2 , where K{x) stands for the com¬ 
plete elliptic integral of the first kind & The numerical 
evaluations of J for the case T = 0 are shown in Fig. [2] 
Considering a 5-nm GaAs coupled quantum well (QW), 
separated by a 4-nm Alo. 33 Gao. 67 As barrier [32) with the 
corresponding effective electron mass m e = 50.061?no 
(to 0 is the free-electron mass) and in-plane effective hole- 
mass of around O.lnio (according with the Luttinger pa¬ 
rameters [ 32 !), we use the exciton binding energy (§2, Slj 
value 6.7 meV to calculate the energy scales in the model. 
Particularly, for the energy scale corresponding to J we 
find for the quasi-2D GaAs/AlGaAs QW structure ge¬ 
ometry J » 0.001796 meV (corresponding to U = 0.0188 
eV) or, in temperature units J ft 20 mK. 


V. BKT TRANSITION 

To understand the nature of the 2D excitonic low- 
temperature phase we describe the excitonic system far 
beyond the mean field description SI- We will use the 
forme of phase-stiffness action in Eci. (l2Ul) to derive the 
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harmonic action in the SCHA schem. Furthermore, using 
the Feynmann-Kleinert (FK) method |34t, j35f , we reduce 
the quantum rotor model given by the action in Ea. fllSI) 
to a simple classical Flamiltonian description. First, let’s 
simplify the form of the phase-stiffness action in Eci. (l20l) 
by transforming the linear trigonometric function inside, 


Here, T = F (r — r') is a new stiffness constant, which 
we will determine by using the variational approach in 
the theory of SCHA. We obtain the system of coupled 
SC equations for the bandwidth renormalization factor 
g b given in Eq. m and the constant T in Eq. (l23l) is (the 
case T = 0 is considered here) 



0 1 2 3 4 5 6 7 


gB = 1 - 


TU 1 — cos (kd) 
4 Jg B N^ y/WTjk ’ 


F( r 


H) = 2g B J • exp 


1 ^ 1 — cos (kd) 

2 

(24) 


For a given U and J, the system of SC equations in 
Eci. (l24l) has many roots of g B . Among those solutions, 
only the largest one must be chosen, because it minimizes 
the effective action in Eq. d^l) . By solving numerically 
the system in Ea. (l24l) . we got the largest root of g B be¬ 
ing equal to unity for all fixed values of the pair ( U , J), 
thus gB = 1 in our model (this is true, if we are interest¬ 
ing of the low-temperature bound-vortex-antivortex type 
configuration, but in general g B differs from unity in 
the unbinding-vortex state of the system). The partition 
function of the system is now 


u/t 


FIG. 3: Critical temperature of the BKT-superfluid phase 
transition in the 2D excitonic system. 

into a quadratic one, and by linearizing the latter one 
with using the standard decoupling procedure, thus, we 
have cos 2 [ p ( yt ) — ^(r'r)] = 2 cos 2 [</?(rr) — </j(r'r)] — 1 ss 
4cos [</?(rr) — <p(r'r)] (cos [95 (rr) — ^(r'r)]} — 1. We get 
the linearized version of the action in Eci. (l20l) Sj \<p\ = 
—J Jq drJ2( T r') cos [v 3 ( rr ) ~ V 5 ( r, ' r )]> where J now, is the 

renormalized effective stiffness-parameter J = 2Jg B . 
The total phase action, including also the kinetic term 
is 

r/3 

+J Y^ / d,T {1 — cos \p (rr ) — p (t't )]} . (22) 
(ry) Jo 

We have added a constant term to the phase stiffness 
action Sj [ 95 ], in order to follow the standard notations 
[ 27 L [28]. Next, we apply the_ self-consistent-harmonic- 
approximation (SCHA) [23, |28[ to replace the action in 
Eq. (1^1) by a trial harmonic action Sscha M 




— J Y, cos [p{v r) — ip(r'r)\ 


(r.r') 


(25) 


Now, we will derive from Ea. (l25l) the effective classical 
Hamiltonian analogue to that of 2D X — Y model. We 
use the FK method for the classical action in Eci. (l25l) 
combined with the SCHA method. Following [34(, we 
decompose the phase variables into static and periodic 
parts p (rr) = p 0 (r) + [<Pn (r) e“" r + c.c] . Then, 

using the FK method, we choose a trial partition func¬ 
tion Z' corresponding to Z in Eci. (l!?5ll . in which the 
potential energy is approximated by a Gaussian com¬ 
ponent for the part of the action with n 7 ^ 0. We 
employ the extremal principle to find Z' and we have 
= JUPV] e where the classical effective 

Hamiltonian TL e s [p\ is given by 

'HeS [<p] = —-7"eff [p\ — 

-7TJ Y ° 2 ( r ) / (D r ') C0S [^( r ) - ^( r ')] - 
(r,r'> 

-2?r J Y f ( r > r ') cos b( r ) - ¥>(0] ■ (26) 

<r,r'> 


Here, the free energy jF e s [p\ is 


•T'eff [<p] = 


1 

0 


E ln 


sinh (j3Un r /2) 
pUtl r /2 


(27) 
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The coefficients / (r, r') stands for the nor¬ 
mal Gaussian smearing factors /( r,r') = 

exp [—a 2 (r) /2 — a 2 (r') /2]/27r. The parameters S2 r 
and a 2 (r) are inter-related by the minimization proce¬ 
dure. We have 


‘ J I / 

f ( r ’ r ') C ° S fa( r ) “ <P( r ')] 


a 2 (r) = 


Ut3VL 2 r 


u/3n r 


coth 




-1 


(28) 


where the prime on the sum in the first of equations in 
Ea. (l28l) is over all r' nearest neighbors of r. Furthermore, 
we have in Eg. (1281) an infinite number of coupled SC 
equations, which do not allow the solution in general. In 
contrast, at the low temperatures we can neglect thermal 
fluctuations of fl r and a (r) and we suppose that fl T = fl 
and a (r) = a, i.e. the uniform solution. Hence, in the 
low-temperature lirnite, we have the classical Hamilto¬ 
nian of the system given by 


Heff M = - Jg cos ( r ) - T ( r 0] , (29) 

( r . r '> 


where g = go (1 — In go) and, for the pair of parameters 
( U , J), go is given by the following equation 


go = 


JU 


T 


) •tanh 1 2 

J 

lu go 

) 

t\ 

1 J 


(30) 


In fact, the classical phase Hamiltonian in Ea. (l28l) is a 
usual X — Y model Hamiltonian in 2D. Considering the 
bound-vortex-antivortex excitation state, the Hamilto¬ 
nian in Eg. (1291) gives rise to a temperature-induced tran¬ 
sition in the long wavelength limit of phase fluctuations. 
The equation for the BKT transition temperature is 

Tbkt = (Tbkt) g- (31) 

Furthermore, Eas. (l30l) and ED- form a system of cou¬ 
pled SC equations for the parameter go and for the BKT 
transition temperature Tbkt- We should take into ac¬ 
count only the higher values of the parameter go and the 
temperatures Tbkt corresponding to it, for all other fixed 
values of the pair (U, J). Using again the fast convergent 
Newton’s method, [291 ] with the relative error of order 
10 -7 , we find the solution (go, Tbkt) of the system. The 
corresponding values of the critical temperature Tbkt 
of the BKT superfluid transition are given in Fig. [3] for 


two different values of the /-band hopping t. For energy 
scale of Tbkt in the case of quasi-2D GaAs/AlGaAs QW 
structure geometry (see the Section [IV] for structure de¬ 
tails) we find Tbkt ~ 12.06 nK at U = 0.0188 eV. For a 
Ino. 08 Gao. 92 As/GaAs coupled QW’s the BKT transition 
temperature is slightly higher and Tbkt ~ 15.14 nK, at 
the same values U = 0.0188 eV. Note, that the recent 
experimental investigations for the probe of superfluid¬ 
ity in 2D Bose systems EH have proved the existance of 
the Landau critical velocity in the temperature range of 
order 90 nK. 


VI. CONCLUSIONS 

We have shown the existence of the BKT phase tran¬ 
sition in 2D EFKM. Using the path integral approach 
and U(l) gauge representation of electron operator, we 
derived the general form of the action of the system. 
Then, we integrated out the fermionic and bosonic de¬ 
grees of freedom, and we got the effective phase action 
(for the bosonic sector) and effective fermionic action (for 
the fermionic sector). Furthermore, the series expansion 
of the effective phase action up to second order in hop¬ 
ping amplitude, gives the form of the parameter J of 
the phase-stiffness between n.n. excitonic pairs. Then, 
in the limit of low-temperatures, we have considered the 
“vortex-antivortex” bound state modes, and we have re¬ 
duced the quantum rotor model to the simple classical 
Hamiltonian description with the temperature dependent 
effective phase-stiffness parameter. Then, we calculated 
the BKT transition critical temperature Tbkt by solv¬ 
ing numerically a system of coupled SC equations for 
the parameters g and Tbkt- Thereby, we have shown 
that the BKT transition in the system of preformed ex¬ 
citonic pairs is directly related with the bound-vortex- 
antivortex type phase-configuration in the system, and 
the excitonic insulator state, with the order-parameter 
A, is a necessary prerequisite for the realization of such 
a transition. In the frame of the same methods we used 
here, it still fundamental to answer the question, whether 
excitonic BEC transition temperature is coinciding with 
the critical temperature Tbkt of the excitonic superfluid 
transition, and for this, the quasi-2D system of excitons 
should be considered and the inter-layer exciton corre¬ 
lations would be properly included. We will discuss on 
that subject in the near future. 
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